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ABSTRACT 

Wc compute electromagnetic wave propagation through the magnetosphere of a mag- 
netar. The magnetosphere is modeled as the QED vacuum and a cold, strongly mag- 
netized plasma. The background field and electromagnetic waves are treated nonpcr- 
turbatively and can be arbitrarily strong. This technique is particularly useful for 
examining non-linear effects in propagating waves. Waves travelling through such a 
medium typically form shocks; on the other hand we focus on the possible existence 
of waves that travel without evolving. Therefore, in order to examine the nonlinear 
effects, we make a travelling wave ansatz and numerically explore the resulting wave 
equations. We discover a class of solutions in a homogeneous plasma which are sta- 
bilized against forming shocks by exciting nonorthogonal components which exhibit 
strong nonlinear behaviour. These waves may be an important part of the energy 
transmission processes near pulsars and magnetars. 



1 INTRODUCTION 

■ The magnetosphere of a magnetar provides a particularly 
, interesting medium for the propagation of electromagnetic 
' waves. Magnetars are characterized by exceptionally large 

magnetic fields that can be several times larger than the 
, quantum critical field strength (Mereghetti 2008). Because 

■ the magnetic fields are so large, the fiuctuations of the vac- 
uum of quantum electrodynamics (QED) influence the prop- 

. agation of light. Speciflcally, the vacuum effects add nonlin- 
' ear terms to the wave equations of light in the presence 
, large magnetic fields. In addition, the magnetosphere of a 

■ magnetar contains a plasma which alters the dispersion re- 
' lationship for light. Because of the unique optical conditions 

in the magnetospheres of magnetars, they provide excellent 
arenas to explore nonlinear vacuum effects arising due to 
quantum electrodynamics. 

The influence of QED vacuum effects from strong mag- 
netic flelds in the vicinities of magnetized stars has previ- 
ously been studied by several authors. The combined QED 
vacuum and plasma medium is discussed in detail in the 
context of neutron stars in Meszaros (1992). Vacuum effects 
have been found to dominate the polarization properties and 
transport of X-rays in the strong magnetic fields near neu- 
tron stars (Meszaros & Ventura 1979, 1978; Meszaros & 
Bonazzola 1981; Meszaros et al. 1980; Galtsov & Nikitina 
1983). Detailed consideration of magnetic vacuum effects 
is therefore critical to an understanding of emissions from 
highly magnetized stars. 

Most studies of waves in systems including plasmas or 
vacuum effects approach the problem perturbatively, which 
limits the applicability of their results. The purpose of the 
present paper is to examine the combined impact of the 



QED vacuum and a magnetized plasma using nonpertur- 
bative methods to fully preserve the nonlinear interaction 
between the fields. Neutron stars may be capable of produc- 
ing very intense electromagnetic waves, comparable to the 
ambient magnetic field. For example, a coupling between 
plasma waves and seismic activity in the crust could pro- 
duce an Alfven wave with a very large amplitude (Blaes 
et al. 1989; Thompson & Duncan 1995). Even if they may 
not be produced directly, electromagnetic waves naturally 
develop through the interactions between Alfven waves. To 
lowest order in the size of the wave, Alfven waves do not suf- 
fer from shock formation (Thompson & Blaes 1998) whereas 
electromagnetic waves do (Heyl & Hernquist 1998, 1999a); 
therefore, how to stabilise the propagation of the the latter 
is the focus of this paper. 

If such a magnetospheric disturbance results in elec- 
tromagnetic waves of sufficiently large amplitude and low 
frequency, then nonperturbative techniques are required to 
characterize the wave. The importance of studying such a 
system nonperturbatively is particularly well illustrated by 
the fact that some nonlinear wave behaviour is fundamen- 
tally nonperturbative, as is generally the case with soli- 
tons (Rajaraman 1982). In order to handle the problem non- 
perturbatively, we choose to study waves whose spacetime 
dependence is described by the parameter S = x — vt, where 
D is a constant speed of propagation through the medium in 
the x-direction. In the study of waves, one normally chooses 
the ansatz e''"*^''"\ However, in this picture, a numeri- 
cal study would typically treat the self-interactions of the 
electromagnetic field by summing the interactions of finitely 
many Fourier modes. So, this ansatz confiicts with our goal 
of studying the nonlinear interactions to all orders. In con- 



trast, a plane wave ansatz given hj S = x ~ vt allows us 
to study a simple wave structure to all orders without any 
reference to individual Fourier modes. 

We model a magnetar atmosphere by including the ef- 
fects of arbitrarily strong electromagnetic fields using a QED 
one- loop effective Lagrangian approach. These effects are 
discussed in section 2.2. Plasma effects are included by as- 
suming free electrons moving under the Lorentz force with- 
out any self-interactions. The model is that of a cold mag- 
netohydrodynamic plasma and is discussed in section 2.3. 
We have also assumed that the medium is homogeneous in 
agreement with the travelling-wave ansatz. Of course the ac- 
tual situation is more complicated with a thermally excited 
plasma (e.g. Gill & Heyl 2009) and inhomogeneities — the 
latter can result in a whole slew of interesting interactions 
between the wave modes (Ifeyl & Shaviv 2000, 2002; Heyl 
et al. 2003; Lai & Ho 2003) that are especially crucial to our 
understanding of the thermal radiation from their surfaces, 
but these are beyond the scope of this paper. 

The formation of electromagnetic shocks is expected 
to be an important phenomenon for electromagnetic waves 
in the magnetized vacuum since electromagnetic waves can 
evolve discontinuities under the influence of nonlinear in- 
teractions (Lutzky & Toll 1959; Zheleznyakov & Fabrikant 
1982; Heyl & Hernquist 1998). Such shocks can form even 
in the presence of a plasma (Heyl & Hernquist 1999b). In 
this study, through our explicit focus on travelling waves, 
we examine an alternate class of solutions to the wave equa- 
tions which do not suffer this fate. Instead, they are stabi- 
lized against the formation of discontinuities by nonlinear 
features. These waves travel as periodic wave trains with- 
out any change to their form, such as wave steepening or 
shock formation. Waves such as these may contribute to the 
formation of pulsar microstructures (Chian & Kennel 1987; 
Jenet et al. 2001). 



2 WAVE EQUATIONS 

2.1 The Maxwell's equations 

The vacuum of QED in the presence of large magnetic fields 
can be described as a non-linear optical medium (Heyl & 
Hernquist 1997b). We also choose to treat the effect of the 
plasma on the waves through source terms pp and Jp; there- 
fore, we begin by considering Maxwell's equations in the 
presence of a medium and plasma sources. In Heaviside- 
Lorentz units with c = 1, Maxwell's equations can be used 
to derive the wave equations 



-V X (V X (D- E)) -f 
^(Vx(B-H)) + Vpp+^ (1) 



-V(V-(B-H)). 
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dt 



(Vx (D-E))-Vx Jf 



(2) 



For clarity, we will avoid making cancellations or dropping 
vanishing terms. We define the vacuum dielectric and inverse 
magnetic permeability tensors as follows (Jackson 1975) 



D, = Eij-Ej, Hi = Bj. 



(3) 



In the next few sections we build a model describing travel- 
ling waves in a magnetar's atmosphere from these equations. 



2.2 Vacuum Dielectric and Inverse Magnetic 
Permeability Tensors 

In this section, we describe our model of the QED vacuum in 
strong background fields in terms of vacuum dielectric and 
inverse magnetic permeability tensors. These are most con- 
veniently described in terms of two Lorentz invariant com- 
binations of the fields. 
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In order to examine the nonlinear effects of the vacuum non- 
perturbatively, we wish to use vacuum dielectric and inverse 
magnetic permeability tensors which are valid to all orders 
in the fields. Analytic expressions for these tensors were de- 
rived by Heyl & Hernquist (1997b) for the case of wrenchless 
fields {K = -(4E ■ B)^ = 0) from the Heisenberg-Euler- 
Weisskopf-Schwinger (Heisenberg & Euler 1936; Weisskopf 
1936; Schwinger 1951) one-loop effective Lagrangian in Heyl 
& Hernquist (1997a) and expressed in terms of a set of an- 
alytic functions. 
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The tensors we need are derived in Heyl & Hernquist 
(1997b), except that we have kept terms up to linear or- 
der in the expansion about K = Q instead of dealing with 
the strictly wrenchless case. Our analysis therefore requires 
that K Bt. 
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the fine-structure constant is a = e^/47r in tliese units where 
we have set ?i = c = 1, and 



J/(Pp)v 



d-Ez. ig the mean electron density 



where Vp 

in the plasma in the reference frame where Vp = 0. It is 
important to note the distinction between the mean plasma 
speed, Vp, and the propagation speed of the wave, v. The 
Lorentz factor 7 = _,^ accounts for a relativistic increase 



in the plasma density since cJ'^x = 7^^ci^x' 

Next, we take a time derivative of the current and ex- 
press this in terms of a three-force acting on the plasma. 



-Bfc V 2 



(14) 



Equations (11) and (12) define our model for the QED vac- 
uum in a strong electromagnetic field. Because we will focus 
on photon energies much lower than the rest-mass energy 
of the electron, we have treated the vacuum as strictly non- 
linear. It is not dispersive. The treatment of the dispersive 
properties of the vacuum would require an effective action 
treatment (e.g. Cangemi et al. 1995a,b) rather than the local 
effective Lagrangian treatment used here. 
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Noting that Jp is a function only of 5, we insert the 
Lorentz force F = e(E + Vp x B) and arrive at an expression 
that can be substituted into the source terms, equations (16) 
and (17), to find 3p{S). 



2.3 Plasma 

To investigate travelling waves, we choose our coordinate 
system so that the x-direction is aligned with the direction 
of propagation. Then, the spacetime dependence of the fields 
and sources is given by a single parameter S = x — vt where 
V is the constant phase velocity in the x-direction of the 
travelling wave. At this point, we are choosing to work in a 
specific Lorentz frame. 

We model the plasma as a free electron plasma which 
enters the wave equation through the source terms Pp and 
Jp. For electromagnetic fields obeying the travelling wave 
ansatz, the sources must also obey the ansatz. Then, pp and 
Jp are functions only of S. We therefore treat them as an 
additional field which is integrated along with the electro- 
magnetic components of the field. 

In order to perform the numerical ODE integration for 
the currents, we wish to find expressions for 
For travelling waves, the continuity equation is 

where we are using index notation to label our explicitly 
cartesian {x, y, z} coordinate system. Repeated indices are 
summed. However, whenever x or z appears as an index, it 
is fixed and does not run from 1 to 3. 

We use equation (15) to rewrite the source terms from 
equation (1) 
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Similarly, the source term from equation (2) is 
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To find an expression for we express the current as 
an integral over the phase-space distribution of the electrons 
and linearize the plasma density. 
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The equations we have given above describe a cold, relativis- 
tic, magnetohydrodynamic plasma. Forces on the plasma 
arising due to pressure gradients and gravity are neglected. 
Moreover, in our simulations, we neglect the forces on the 
plasma due to the magnetic field. This approximation is suit- 
able if the plasma in question is a pair plasma, or for wave 
frequencies much less than the cyclotron frequency. 

As we will show in section 4 (see figure 6), the field con- 
figurations generated in our simulation vary over timescales 
similar to the inverse of the plasma frequency; the latter is 
about 9 orders of magnitude smaller than the cyclotron fre- 
quency. This observation allows us to justify some aggressive 
assumptions regarding the plasma response. As mentioned 
above, we may neglect forces on the plasma from the mag- 
netic field, and quantum effects are expected to be small far 
away from the cyclotron resonance (Meszaros 1992). 

We also neglect thermal effects since the influence of the 
electromagnetic fields will dominate over thermal motion. 
For strong background magnetic fields such that — ^ kT, 
the electrons will be confined to the lowest Landau level, 
restricting thermal motion perpendicular to the background 
magnetic field. As we will elaborate on in section 4, the 
greatest nonlinear effects occur for waves with an electric 
field component oriented along the background magnetic 
field. Because we are interested in waves with amplitudes 
comparable to Bk, thermal motion is negligible relative to 
the dynamics induced by the wave. We are therefore justi- 
fied in neglecting thermal effects in every direction for the 
cases of greatest interest. 

If one combines Eq. (19) with Eq. (1), one sees that 
any nonlinearity in this treatment must originate with the 
dielectric and permeability tensors — any non-linearity that 
may originate from the plasma itself has been neglected (c.f. 
Kozlov et al. 1979; Cattaert et al. 2005). The plasma is mod- 
elled as strictly dispersive. In section 2.5, we confirm that 
this method of describing the plasma is consistent with stan- 
dard accounts in the weak-field, small- wave limit. 



2.4 Travelling Wave ODEs 

The wave equations for travelling waves, are found by com- 
bining Maxwell's equations (1) and (2) with the continuity 
equation, (15). We also make the plane-wave approximation, 
and assume that the fields and sources are described by the 
parameter 



S^x-vt. (21) 
The equations governing travelling wave propagation are 
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The auxiliary vectors and are related to the electric 
and magnetic fields. 
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Equations (22) through (25) define a set of coupled ordinary 
differential equations that can be integrated to solve for the 
travelling electric and magnetic fields. Solving these equa- 
tions requires that we have at hand the vacuum dielectric 
and inverse magnetic permeability tensors as well as an ex- 
pression for These were discussed in sections 2.2 and 
2.3 respectively. 



2.5 Weak field, Small Wave limit 

In this section, we would like to demonstrate that our equa- 
tions reduce to standard expressions in the case of small 
background fields and small electromagnetic waves. To make 
this comparison, it is also prudent to assume waves have a 
spacetime dependence like e*'"'"''"' instead of a: — vt. By 
making this change, we can compare our other assumptions 
with those made in standard textbook accounts directly. 

Under the assumption that the fields and currents have 
the standard plane-wave spacetime dependence, we can 
make the replacements ^ — >• —iuj and V — > ik. We can 
then write a second expression for — 



dt ■ 



dJ 



Setting equations (26) and (19) equal, we get 
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where we have used the nonrelativistic approximation, which 
is appropriate for small waves. We will now specialize to a 
background magnetic field pointing in the z-direction. We 
can then write this background field in terms of the cyclotron 



frequency ojp = Then, we can write equation (27) as 
[-LUc£*'" - iujS''^ 4 = ujIe' (28) 
We can solve this for Jp by writing a matrix equation 




inverting the matrix gives 
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Finally, we would like to use equation (30) to write the 
right hand side of equation (1) in the small wave limit in a 
manner we can interpret as a dielectric tensor. 

Ignoring (for now) the contributions from the vacuum, 
and assuming an approximately homogeneous plasma den- 
sity, equation (1) simplifies to 



dt 

— iojJr: 



(31) 



If our macroscopic field is to obey V^D — = 0, we 

can insert equation (30) into (31) to obtain the following 
expression for the dielectric tensor due to plasma effects: 
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This expression is in agreement with the cold plasma dielec- 
tric tensor given in Meszaros (1992). As noted above, our 
analysis neglects the off-diagonal (Hall) terms, as is appro- 
priate for pair plasmas or waves with frequencies much less 
than the cyclotron frequency. 

Because the vacuum effects are added explicitly in the 
form of dielectric and magnetic permeability tensors, we only 
need to confirm that the weak field limits of our expressions 
agree with the standard results. This confirmation is done 
explicitly in Heyl & Hernquist (1997b). 

In the weak field limit, the tensors given by equations 
(11) and (12) are 

^ " [2{E^ - B'^)Sij + TBiBj 
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In the case of a weak background magnetic field pointing in 
the z-direction, these become 
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The density n determines the plasma 
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Figure 1. A comparison between tlie z components of the electric 
(solid) and magnetic fields (dashed) showing the nonorthogonal 
stabilizing wave. On larger scales, the component is seen to 
have a periodic envelope structure as in figure 2 
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Figure 2. The wave component as shown in figure 1 has an 
envelope structure when viewed at larger scales. In this case, the 
Ez wave component has a nearly constant amplitude of O.IB],. 



Again, this result agrees with the vacuum tensors given in 
Meszaros (1992). 

In this hmit, we may simply add together the contri- 
butions to the dielectric tensor from the plasma and the 
vacuum according to 

e^,^5.,+{e^' + -hj) (38) 

with given entirely by the vacuum contribution. 

We have thus recovered the standard result for a 
medium consisting of a plasma and the QED vacuum in 
the weak field, small wave limit. 



3 SOLUTION PROCEDURE 

In total, there are 15 coupled non-linear ODEs which must 
be integrated to produce a solution. Equations (22) and (23) 
define the electric and magnetic fields as functions of 5*. In 
addition, we must simultaneously integrate equation (20) 
which gives the plasma current as a function of S. Each of 
these equations has three spatial components. Initial condi- 
tions are given for each of the six field components, and the 
six derivatives of the field with respect to S. The equations 
describing the electromagnetic fields do not depend on the 
initial values of the current, but in principle these are the 
three remaining initial conditions. 

The ODEs are solved using a variable stepsize Runge- 
Kutta method. In order to translate between the E-and-B 



picture and the ?/)-and-x picture, equations (24) and (25) 
must be solved numerically at each time step, including the 
first step when the initial conditions are given. 

Tables of numerical values of the functions defined by 
equations (6) to (8), as well as their derivatives, were com- 
puted in advance and these were used to interpolate the val- 
ues needed in the simulation using a standard cubic spline 
interpolation algorithm. In producing these tables, expres- 
sions for the weak and strong field limits were used in the 
appropriate regimes as this reduced the numerical errors. 

Aside from the initial conditions for the fields and 
derivatives, there is one parameter in the model which 
must be selected. The density of the plasma, given by n 
in equation(18) is chosen to be 

n = lO^^cm"^ (39) 

This value corresponds to the Goldreich- Julian density for a 
star of period P ~ 1 s, P ~ 10^^". The uniform background 
field is taken to equal the quantum critical field strength 

2 

Bfc = — = 4.413 X lO'^G. (40) 



4 RESULTS 

This study focuses on the case of waves propagating trans- 
verse to a large background magnetic field. We have already 
chosen the direction of propagation to be the x-direction 
through our definition of S in Eq. (21). We now choose the 
background magnetic field to point in the z-direction. In this 
situation, the largest non-linear effects occur when there is 
a large amplitude wave in the z component of the electric 
field. This is quite natural. The values of both the dielectric 
tensor of the plasma (Eq. 32) and the weak-field limit of 



the vacuum dielectric tensor (Eq.35) differ most from unity 
for this component; for strong fields the index of refraction 
(as well as its derivative with respect to the field strength) 
is largest for vacuum propagation in this mode (Heyl & 
Hernquist 1997b). Furthermore, the dominant three-photon 
interaction couples photons with the electric field pointed 
along the global magnetic field direction with photons whose 
magnetic field points along this direction. The three-point 
interaction for photons whose electric field is perpendicular 
to the magnetic field with parallel photons vanishes by the 
CP— invariance of QED (Adler 1971). Therefore, we focus 
on initial (5* = 0) conditions in which the dominant compo- 
nent of the electric field points along in the z-direction and 
that of the magnetic field along the y-direction. 

In the classical vacuum, these initial conditions corre- 
spond to transverse, linearly polarized sine-wave solutions 
that travel at the speed of light. However, when the wave 
amplitudes are large, and there is a strongly magnetized 
plasma, we find that there is a deviation from normal trans- 
verse electromagnetic waves. In particular, in order to re- 
main stable, a wave with large and By field components 
must also excite waves in the Ey and fields. These sta- 
bilizing wave components exhibit strong non-linear charac- 
teristics (see figures 1 and 2). The symmetries of the wave 
equation require a close correspondence between the Ey and 
Bz waveforms as well as between the Bz and By waveforms. 
For simplicity, only one of each is plotted in the figures. The 
field strengths are given in units of the quantum critical field 
strength, Bk- 

As is apparent from Fig. 1 the dominant electric field 
along the direction of the external magnetic field is essen- 
tially sinusoidal. Subsequent figures will show that there is a 
small harmonic component. The waveform for the dominant 
magnetic field component is similar. On the other hand, the 
magnetic field along the direction of the electric field (the 
non-orthogonal component) is smaller by nearly four orders 
of magnitude and obviously exhibits higher harmonics. In 
particular if one expands the scale of interest (Fig. 2), the 
magnetic field exhibits beating between two nearby frequen- 
cies with similar power. 

In order to examine the harmonic content of the wave- 
forms, we perform fast Fourier transforms (FFTs) on the 
signals produced in the simulations. We present the results 
in terms of power spectra normalized by the square am- 
plitudes of the electric field of the waves. In these plots, 
the horizontal axis is normalized by the frequency with the 
greatest power in the electric field, so that harmonics can be 
easily identified. 

Fig. 3 depicts the power spectra of the electric and mag- 
netic fields along the z-direction for the wave depicted in 
Figs. 1 and 2. The conclusions gathered from an examina- 
tion of the waveforms are born out by the power spectra. In 
particular the electric field is a pure sinusoidal variation to 
about one part in ten thousand - the power spectrum of a 
pure sinusoid is given by the dashed curved. The duration 
of the simulation is not an integral multiple of the period 
of the sinusoid, resulting in a broad power spectrum even 
for a pure sinusoid. The power spectrum of the magnetic 
field follows the expectations gleaned from the waveforms. 
In particular the fundamental and the first harmonic are 
dominant, with the first harmonic having about one-third 
the power of the fundamental. If one focusses on the fun- 
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Figure 3. The solid curve depicts power spectra of the electric 
(upper) and magnetic fields (lower) along the z-direction for the 
solution depicted in Figs. 1 and 2. The inset focusses on the fun- 
damental and the first harmonic. The dashed curve follows the 
power spectrum of a single sinusoid for the electric field and three 
sinusoids for the magnetic field. Near the peaks the dashed curve 
is essentially indistinguishable from the solid one. 



damental, one sees that two frequencies are involved. The 
envelope structure is produced by a beating between the 
fundamental and a slightly lower frequency with a similar 
amount of power as the first harmonic. Over the course of 
the simulation the envelope exhibits two apparent oscilla- 
tions; the lower frequency differs by two frequency bins, so 
it is resolved separately from the fundamental as shown in 
the inset. 

As the amplitude of the electric field increases the non- 
linear and non-orthogonal features of the travelling wave in- 
crease. Fig. 4 shows that the strength of the non-orthogonal 
magnetic field increases as the square of the electric field, a 
hallmark of the non-linear interaction between the fields. For 
the strongest waves studied with E^ ~ 0.2Bk (the rightmost 
point in the figure), the magnetic field, B^, is about 10""^ Bk 
nearly one-percent of the electric field. The amplitude of the 
non-orthogonal magnetic field is given by 



Bz = O.OOSSfc ( ^ 
, i>fc 



(41) 



for Bo = Bk- The coefficient is coincident ally very close to 
three-quarters of the value of the fine-structure constant. It 
increases with the strength of the background field. 

For the strongest waves even the non-orthogonal mag- 
netic field is strong, so it can generate non-linearities in the 
electric field. Although the strongest effect is around the 
fundamental, it is completely swamped by the fundamental 
of the electric field. On the other hand, the magnetic field 
drives a first and second harmonic in the electric field as seen 
in Fig. 5. The strength of these harmonics is approximately 
given by the formula in Eq. (41) or equivalently Fig. 4 if one 
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Figure 4. Amplitude of tiic Bz component plotted against the 
amplitude of the component for a Bq = Bf^ background field. 
The line is the best-fitting power-law relation. The slope is con- 
sistent with a scaling exponent equal to two. 



substitutes the value of Bz for Ez and uses result for Ez. 
This is essentially a six-order correction from the effective 
Lagrangian. Because we have used the complete Lagrangian 
rather than a term-by-term expansion, all of the corrections 
up to sixth order (and further) are automatically included 
in the calculation. 

Figure 6 demonstrates how the frequency of the solu- 
tions varies with the speed of propagation. Because there is 
no mode information stored directly in our numerical solu- 
tions, we take the frequency to be rate that local minima in 
the electric field pass a fixed observer. In general, we find 
that the frequency of traveling waves increases as the phase 
velocity approaches the speed of hght, very closely following 
the formula 



1.2 - 1 



(42) 



This formula also results from an analysis of the dielectric 
tensor, Eq. (32) . The vacuum makes a small contribution to 
the wave velocity in this regime. 



5 CONCLUSION 

We have discussed techniques for computing electromagnetic 
waves in a strongly magnetized plasma which nonperturba- 
tively account for the field interactions arising from QED 
vacuum effects. We applied these methods to the case of 
travelling waves, which have a spacetime dependence given 
by the parameter S — x — vt. Travelling waves can be de- 
scribed without any decomposition into Fourier modes and 
this is ideal for exploring the nonperturbative aspects of 
waves. 

The main result from this analysis is the observation 
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Figure 5. This power spectrum demonstrates the development 
of nonlinear effects in the z-component of the electric field as the 
amplitude of the wave is increased in a Bq = Bj, background 
field. From top to bottom the curves follow the solutions whose 
amplitude of Ez equals 0.01,0.08 and 0.16-Bfe. The inset focusses 
on the first and second harmonics. 
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Figure 6. The frequency (solid curve) and wavenumber (dashed 
curve) of the travelling wave as a function of its phase velocity. 



that electromagnetic waves in a strongly magnetized plasma 
can self-stabilize by exciting additional nonorthogonal wave 
components. In the cases studied above, a large amplitude 
excitation of the electromagnetic field, for example, from the 
coupling between Alfven waves to starquakes, can induce 
nonlinear waves which are stabilized against the formation 



of shocks. The result is a periodic wave train with distinctly 
nonlinear characteristics. Such structures may play a role in 
forming pulsar microstructures. 

This result demonstrates that shock formation is not 
a necessary outcome for waves in a critically magnetized 
plasma. It is possible that nonlinear features of a wave can 
stabilize it against shock formation. The shock-wave solu- 
tions generally decrease in magnitude, dissipating energy 
along the way unless they travel into a low-field region, so 
the self-stabilizing nonlinear waves are the only ones that 
keep their shape and energy content intact as they propa- 
gate. It is not yet clear what set of conditions will determine 
if a particular wave will self-stabilize or collapse to form a 
shock or how plasma inhomogenaities will affect the propa- 
gation of a travelling wave train. These are issues which can 
be clarified in future work. 
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